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Abstract 

We study the time evolution of classical spin systems with purely relaxational 
dynamics, quenched from T ^ Tc to the critical point, in the semi-infinite 
geometry. Shortly after the quench, like in the bulk, a nonequilibrium regime 
governed by universal power laws is also found near the surface. We show for 
'ordinary' and 'special' transitions that the corresponding critical exponents 
differ from their bulk values, but can be expressed via scaling relations in terms 
of known bulk and surface exponents. To corroborate our scaling analysis, we 
present perturbative (e-expansion) and Monte Carlo results. 
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The dynamics of spin systems after quenches from a high-temperature initial state to 
or below a critical point has attracted a great deal of interest during recent years 0]. In 
particular, for quenches from T ^ to the critical temperature Janssen et al. found that 
in Ising-like systems with purely relaxational dynamics (model A) a previously unknown 
regime shortly after the quench with unexpected features exists. In its most pronounced 
form the so-called universal short-time behavior (USTB) appears in the time dependence 
of the order parameter itself. Suppose the initial state carries some small magnetization 
mo. Then, as Janssen et al. showed, after the quench the magnetization will grow for a 
macroscopic time span, governed by the universal power law 

m(t)~mot^ with 6 = {xq - l3/iy)/C, (1) 

where Xq is the scaling dimension of the initial field (whereas the static exponents have 
their usual meaning, and ( denotes the dynamic exponent). Thus, technically spoken, USTB 
is caused by the influence of the scaling behavior of the initial field on (in macroscopic units) 
early stages of the relaxational process [0. The corresponding operator dimension xq is a 
genuine new exponent. Especially for the n- vector model in 2 < (i < 4 it turns out that 
Xq > jS/u, and thus the magnetization really increases initially, up to a characteristic time 
to ~ tuq'^^'"'. After that it crosses over to the well-known algebraic decay m t '^'^^ [PIH- 
Recently, the scenario developed in Ref. [0] was directly confirmed by means of Monte Carlo 
simulation |]^, and, based on USTB, a procedure to determine static critical exponents was 
suggested 0. 

Another quantity where the short-time exponent 9 enters is the autocorrelation function 
A{t) = (0(t) 0(0)) ~ t"'^/''. It was shown by Janssen that \ = d ~ (6, and it is this 



quantity where USTB was actually first seen in a numerical simulation . 

Similar in many respects to initial conditions for time-dependent processes are boundary 
conditions at surfaces. It is well known that critical behavior is modified near surfaces; each 
bulk universality class splits up in three surface universality classes, depending on whether 
at bulk criticality the surface is disordered, ordered, or also critical. The corresponding 
universality classes are called ordinary, extraordinary, and special transition, respectively 
|ll|]. In this Letter we focus attention on the influence of 'ordinary' and 'special' surfaces on 



the time-dependence of m and A. We show that close to the surface the USTB is modified 
depending on the surface universality class. Further, we show that the values of the 'new' 
exponents that govern the short-time dynamics near the surface can be expressed in terms 
of scaling relations between other known critical exponents. In a recent preprint, Majumdar 
and Sengupta [|1^ correctly determined the surface behavior of the autocorrelation function 



for the n-vector model, but the underlying scaling relation was not recognized. 

In the framework of field theory the modification of critical behavior near surfaces is the 
result of additional primitive uv-divergences, which are located in the surface These 



divergences need additional Z-factors, which, in turn, give rise to independent surface expo- 
nents. As shown by Janssen et al. 0, the situation is largely analogous at the time 'surface' 
t = 0. That in general xo ^ (3/v is the result of an additional logarithmic divergence at t = 
in the upper critical dimension d* = 4. What happens, if we have both initial and boundary 
conditions? It is well known that the degree of surface divergences is reduced at least by one 



compared with bulk divergences Since the bulk divergence at t = is just logarithmic, 



it is immediately clear that no additional primitive divergence located at z = t = exists, 
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and, thus, the USTB near surfaces must be governed by known bulk, bulk short-time, and 
surface exponents. 

To be more specific, let us consider the t- and z-dependent magnetization in the 
semi-infinite system. At the bulk critical point it follows from dimensional analysis and 
renormalization-group (RG) invariance that asymptotically 

m{z, t, mo) ~ r^/^^ y^{z/t^^^, niQ t^^/^ . (2) 

As usual all quantities are made dimensionless with the help of an appropriate power of 
the renormalization mass fi, and we set = 1 afterwards. Because of different scaling 
dimensions, mo is not the initial value of m. Instead, as close as possible to the naive 
expectation, one finds that m ~ mo for t -C to [0 ■ follows that J-'{x, ?/) ~ ?/ for x — cx) 
and 1/ — > 0, and we obtain the known bulk short-time behavior of Eqn. (|I]). 

When do we expect modifications due to the surface? The surface will influence the 
dynamics of a spin when the nonequilibrium (growing) ation length ^(t) ~ t^^'^ has become 
larger than the spin's distance z from the surface. This is analogous to the static situation, 
where surface critical behavior is observed if the equilibrium correlation length is larger than 
z 



TT[| . Now, ii z <^ ^(t), the short- distance expansion (SDE) [p^p 



m{z, t, mo) ~ ^(^i-^)/'^ mi(t, mo) (3) 



should hold, where Pi/i^ is the scaling dimension of the surface field [|T^ and mo ho- 
mogeneous. Moreover, exploiting as before the RG- invariance leads to mi(t,mo) ~ 
f-l^i/'^zjr^l^^^^xo/zy Demanding again that m ~ mo for t <^ to, it follows that J-'i{y) ~ y for 
y ^ 0. As a result, the magnetization close to the surface {z <^ C,{t)) takes the form 

m(z, t, mo) ~ mo z^^'-^^^" t^' , (4a) 

with 



01 = {xo - Pi/iy)/C = 0+{P- (3r)/K ■ (4b) 

In other words, the short-time exponent near the surface is determined by the difference 
between the scaling dimensions of the initial field and the surface field. Eqn. ( ^a]) holds if 
t <^ to o.iT'd z <^ C,{t). If, on the other hand, z > ^(to) only bulk short-time behavior will be 
seen and eventually a crossover to the nonlinear surface behavior m ~ 

Next we set mo = and consider the autocorrelation function A{z,t) := (0(r, t)(/)(r, 0)). 
Because of the Dirichlet initial condition the asymptotic result for A vanishes identically, 
and, as shown by Janssen one has to consider the leading correction in an expansion in 
the inverse of the initial temperature Tq^: 

A{z,t) = r,\<Pir,t)Mr)). (5) 

00 is the initial response field with bulk scaling dimension d — Xq 0. 

For z <^ C,{t) both operators in (^ can be replaced by their SDE. The SDE of the order- 
parameter field can be read from Eqn. (|^). The SDE of 0o is obtained as follows: It is 
known from the equilibrium dynamics [|T^,|n|] that the response field needs an additional 

1 /2 

renormalization factor Zi at the surface, the same Z-f actor that is required to renormalize 
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the order-parameter field at the surface. Utihzing our ear her resuh that there are no new 

1 /2 

primitive divergencies at 2; = t = 0, we conclude that an additional also suffices to 
renormalize 0o at the surface. As a consequence, 0o has the same SDE as 0, and the 
autocorrelation function takes the form 

(6a) 

with 

Ai = rf-C^ + 2(/5i-/5)/z/. (6b) 

Eqs. (§) and (||) are the essential results of this Letter. Numerical values for Qi and 
Ai, where we have used the literature values of the exponents involved, are given in Table 1 
(together with other results discussed below). 

In order to confirm Eqs. (^) and (^, we present a first-order e-expansion (one-loop) for 
the equation of motion of the n- vector model A . According to the above discussion, m is 
small for t -C to, and so we may consider a linear approximation (in m). It takes the form 

m + rrm-rm" + — (n + 2)C(2,t)m = 0, (7) 

6 

where C(z,t) is the one-loop (tadpole) contribution. According to the standard procedure 
||TT| C{z^ t) can be separated in a bulk and a surface part, where the latter is determined by 
the bulk propagator at the image point: 

C{z,t) =Cb{0,t)±Cb{2z,t). (8a) 

In this equation stands for special and the '— ' for ordinary transition, and 

[2'kY J 

The first term in (Pa|), the bulk contribution, contains the usual uv- divergence. Without 
going into technical details, we state that this divergence (after dimensional regularization) 
can be absorbed in the bare temperature by minimal subtraction as in the static calculation. 
The second term on the right-hand side of (^) is the surface contribution, which is finite 
for any z > 0. So after renormalization we obtain a finite equation, and at the RG fixed 
point g* = 487r^e/(n -|- 8) and r = it takes the form 

m-m" + 9 (^-i ± ^e-"'/*^ m = 0, (9) 



where 6 = e{n + 2)/ (An + 32) is the one-loop result for the bulk short-time exponent 
and with t = 2rt. 

It is straightforward to show that for uniform initial magnetization mo Eqn. (P) has 
solutions of the form m(t, z) = mo U{t) V{z'^/t). Let us first consider the bulk limit z ^ ^{t) 
with ^ ~ t^/^ = 2 in one- loop). By straightforward analysis it follows that U ^ F and 
V const., and thus m . This is the bulk short-time behavior. 

What happens if we move closer to the surface? For z <^ ^(t) the function V{x) behaves 
as a power of a;, which now depends on the universality class of the surface. For the ordinary 
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transition one finds V ~ This means that m ~ uiq ^1-2^^-1/2+29 ^^^^ Qord _ _]^^2 + 

e (n + 2)/(2n + 16). For the special transition the result is ~ and m ~ mQZ~'^^P^ 
such that = e{n + 2)/{2n + 16). The one- loop values of 6°^'^ and 61^ are consistent with 
(^^1). That in this simple approximation the surface exponents can effectively be expressed 
with the help of the bulk 6 is an artifact of the one-loop approximation. 

We do not repeat the calculation of Ref. ||T^ for the autocorrelation function here. The 
one-loop resuhs ^ = 6 - e (5n + 22)/(2n + 16) and Af = 4 - e (5ra + 22) /{2n + 16) derived 
in |12| are both consistent with the scaling relation ([Hbf ). The figures for n = e = 1 are given 
in the Table. 

From the one- loop results we can immediately read off the exact exponents for n —>■ (yo. 
They are given by = {3-d)/2, 6'/ = (4-c^)/2, Af'^ = (5rf-8)/2, and Af = (5rf-12)/2. 
Numerical values for d = 3 are displayed in the Table. 

To summarize these results, at the 'special' surface in 2 < < 4 the order grows always 
faster than in the bulk, and the autocorrelations decay slower. For 'ordinary' surfaces, the 
USTB for the Ising universality class is drastically different from bulk behavior; the order 
starts to decay from the beginning. On the other hand, for n — 00 we find 6"^'^ > for 
d < 3. Autocorrelations generally decay faster at 'ordinary' surfaces. 
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Table 1: Dynamic surface exponents Oi and Xi/C for the d=3 Ising model (n = 1) for 'ordinary' 
and 'special' surfaces. For comparison some bulk results are also displayed. The one-loop approx- 
imation has been extrapolated to e = 1. In section 'MC our Monte Carlo estimates are displayed. 
In 'Lit.' results obtained with ( ^ ) and ( [6b| ) in combination with the best known literature values 
for the exponents involved are shown. The latter are taken from Refs. |1L20,21|. In the last section 
the exact exponents of the n ^ 00 limit in d = 3 are displayed. 



Finally, we present a Monte Carlo simulation of the d = 3 Ising model. We have used the 
heat-bath algorithm for nonconserved spins and the usual time-dependent interpretation of 
our data |1^. The simulations were carried out on a lattice of size 20^ x40, where we have im- 
plemented the nontrivial boundary conditions on the surfaces perpendicular to the long axis 
and periodic boundary conditions in the remaining directions. The high-temperature initial 
configuration was generated as in Ref. []^ by randomly selecting and flipping 1/2 — mo/2 of 
the spins in an initially completely (in -|— direction) ordered configuration. To determine the 
exponent 6 from the order-parameter profiles, we started from an initial state with nonzero 
mo; for the data displayed in Fig. 1 we chose uiq = 0.04. For system size tending to infinity, 
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mo would scale to zero [0. For A(t) a disordered initial state (with mo = 0) was prepared. 
All simulations have been carried out at the bulk critical value J/fc^T^ = 0.2216 [21|. The 



Monte Carlo results for 9i and Ai were obtained by averaging over 200,000 independent 
histories. To estimate the errors of the exponents, we have divided the data into 20 runs. 
To simulate the special transition, the result of Ruge et al. for the critical surface cou- 
pling Ji/J = 1.5004 was adopted. In the case of the ordinary transition we chose a value 
somewhat lower than the bulk coupling ]TP|. 



Fig 1.: Monte-Carlo results from 
20^ X 40 lattice for the order param- 
eter for special (solid curves) and 
ordinary (dashed curves) transition. 
The distances from the surface are 
z = 0, 1, 2, 3, 5, 10, 20. The behav- 
ior at the surface {z = 0) is given by 
the uppermost curve in the case of 
special and by the lowest curve in 
the case of ordinary transition. 



Results for the order parameter for special and ordinary transition are shown in Fig. 1 
for different values of z. As expected from the above discussion, all curves for points off 
the surface start with bulk USTB. Close to the surface the profiles show a crossover to 
the respective surface USTB. Far off the surface this crossover does not occur within the 
period of initial increase. For later times all curves show the 'linear' bulk behavior P,^; the 
nonlinear decay could only be seen in larger systems The values for 6i, determined from 
our data, are displayed in Table 1. Taking into account the errors, they are consistent with 
the expectation from (^). The situation is very similar for the autocorrelation function. 
The results for Ai are also displayed in the Table. 

To conclude, we have investigated the critical dynamics, especially the short-time behav- 
ior, of the n-vector model A near 'ordinary' and 'special' surfaces. By general scaling analysis 
combined with a short-distance expansion, we have related the surface short-time exponents 
to other known exponents. The corresponding scaling relations are given in Eqs. (]^^ and 
(IBBI). In a one- loop calculation these scaling relation were verified, and numerical estimates 
for the surface exponents were determined by Monte Carlo simulation. We think that with 
the help of universal short-time behavior, along the hues suggested recently by Li et al. ^ 
for the bulk, it should also be possible to determine static surface exponents more precisely. 
Further, we believe that short-time dynamics near surfaces could also be an attractive field 
for experimental studies. 
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